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Abstract 

Observational manifestations of some models of modified gravity, which have 
been suggested to explain the accelerated cosmological expansion, are analyzed for 
gravitating systems with time dependent mass density. It is shown that if the mass 
density rises with time, the system evolves to the singular state with infinite curvature 
scalar. The corresponding characteristic time is typically much shorter than the 
cosmological time. 
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Contemporary astronomical data strongly indicate that at the present epoch the uni- 
verse expands with acceleration. A possible way to explain this accelerated expansion is 
to assume that there is a new component in the cosmological energy density, the so called 
dark energy. The latter can be either a small vacuum energy, which is identical to the cos- 
mological constant, or the energy density associated with an unknown, presumably scalar 
field, which slowly varies in the course of the cosmological evolution. 

A competing possibility to create cosmological acceleration is to modify gravity itself 
introducing additional terms into the usual action of General Relativity [Ij; for recent 
reviews see [2,3j. To this purpose the models with the following action were considered: 

S = ^l d'x^gf{R) + Sm. (1) 

where mpi = 1.22 ■ 10^^ GeV is the Planck mass, R is the scalar curvature, and Sm is the 
matter action. In the usual Einstein gravity function f{R) has the form f{R) = R, in the 
modified gravity f{R) acquires an additional term: 

f{R) = R + F{R), (2) 

which changes gravity at large distances and is responsible for cosmological acceleration. 
In the pioneering papers [1] function F{R) at small R behaves as: 

F{R) = -^, (3) 

where is a small parameter with dimension of mass. However, as it was shown in ref. 
such a choice of F{R) leads to a strong exponential instability near massive objects and 
so the usual gravitational fields would be drastically distorted. An attempt to cure this 
ill-behavior by adding to the action gR^-teim [5] was only partially successful. It could 
terminate the instability with reasonably small coefficient g for sufficiently dense objects 
with p > 1 g/cm^, while for the objects with smaller mass density the coefficient g would 
be too large and incompatible with the existing bound on the i?^-gravity. We will discuss 
this problem in a more detailed paper which is under preparation. 

A choice of F{R), which leads to an accelerated cosmological expansion and is devoid 
of the above mentioned instability and of some other problems was suggested in several 
papers [nHH]- In the present work we examine a very interesting model of modified gravity 
with F{R) function suggested in ref. [B] : 



F{R) = \Ro 



1 + — 

Rn 



(4) 



Here constant A is chosen to be positive to produce an accelerated cosmological expansion, 
n is a positive integer, and Rq is a constant with dimension of the curvature scalar. The 
latter is assumed to be of the order of the present day average curvature of the universe, 
i.e. Rq ~ l/tfj, where tu ^ 4 ■ 10^'' sec is the universe age. 
The corresponding equations of motion have the form 

1 O T-'(™) 

(1 + F') R^, --{R + F) g^, + {g^uD^D^ - D^D,) F' = -^^-f- , (5) 
Z rripi 
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where F' = dF/dR, is the covariant derivative, and T^™ is the energy-momentum 
tensor of matter. 

By taking trace over ^ and v in eq. ([5]) we obtain the equation of motion which contains 
only the curvature scalar R and the trace of the energy-momentum tensor of matter: 

SD^F' - R + RF' -2F = T, (6) 

where T = SnTI^ /m1,i. Note that our sign convention is different from that of paper [B] 
and is the same as in ref. [1]. 

Cosmology with gravitational action (jlj), as well as some other cosmological scenarios 
with modified gravity were critically analyzed in recent paper [2]. It was shown that, 
taken literally, the models suffer from several serious problems. Though the instability 
of ref. [1] was eliminated, still there remain some other types of singular behavior. In 
particular, there exists the past singularity, when i? — )■ oo at some finite time in the past. 
It was argued that the problem can be solved by an addition to the action of i?^-term with 
sufficiently small coefficient allowed by the present observational data. 

The singularity similar to that considered in the present work was first noticed in 
ref. [in] in the case of cosmological evolution back to the past. In a sense the future singu- 
larity considered here is the time reversal of the past singularity on the quoted paper. So 
mathematically both singularities are quite similar, despite of some difference due to effects 
of the universe contraction (when one goes backward in time). The Hubble anti-friction 
favors the approach to the singularity in the contracting universe. However, despite math- 
ematical similarities there is an important difference between the two systems. According 
the ref. [TU], the singularity may be avoided with a certain range of initial conditions. In 
our case singularity emerges for any initial conditions. 

In ref. [Ill|T2] it was argued that infinite R singularity could arise in the future, unless 
the initial conditions for R are not fine-tuned. This is similar to the cosmological situation 
of ref. [To]. The systems considered in these works are different from that discussed here 
and the singularity of the quoted papers appears only for certain initial conditions, while 
in our case, as we have already mentioned above, the singularity arises for an arbitrary 
initial state. All the singularities can be eliminated by an addition of i?^-term to the action 
and we study the effects of this term below. 

In what follows we will consider a different physical situation than those discussed in 
the above mentioned references. Namely we study behavior of astronomical objects with 
mass density which rises with time and show that curvature, R, reaches infinitely large 
value during the time interval which is very short in comparison with the cosmological 
time scale. This singularity cannot be eliminated by fine-tuning of the initial conditions. 
An addition of i?^-term could prevent from the singular behavior but at expense of quite 
large values of n which may be at odds with the standard cosmological evolution. 

We study objects with mass density which is much larger than the cosmological one, 
Pm ^ Pc- The cosmological energy density at the present time is pc ~ 10~^^g/cm'^, while 
matter density of, say, a dust cloud in a galaxy could be about pm ~ lO^^'^g/cm^. Since the 
magnitude of the curvature scalar is proportional to the mass density of a nonrelativistic 
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system, we find Rq. In this limit: 



1-1^ 
R 



2n 



(7) 



Let us start from tlie initial state in which modified gravity around or inside some 
massive objects is not much different from the usual Einstein (Newtonian) gravity and 
correspondingly R ^ —T, as can be seen from the normal Einstein equations. 

We analyze temporary evolution of solutions of eq. (E]) for the gravitational field of 
some massive object with time varying density. We assume that the gravitational field 
of this object is weak, as is usually the case. Correspondingly the background metric is 
approximately flat and the covariant derivatives can be replaced by the usual flat space 
ones. Hence: 



D^F' = {dt - A)F' = F'\dt - A)R + F"'[{RY - {VRY] , 

where A is the usual Laplacian, and V is the gradient. 

Substituting expression ([7]) for F{R) at large R into eq. IQ, we obtain: 



R^ 



R 



2n 



3n(2n + 1) V^i 



'2n 




(n + i; 



(9) 



However, due to the presence of the nonlinear terms containing derivatives, this equation 
is difficult to analyze and we instead study the equation for F\R) and express R through 
F' using: 



F' = -2nX 



'r 



2n+l 



(10) 



Notice that infinite R corresponds to F' = and if F' reaches zero, it would mean that R 
becomes infinitly large. 

Let us introduce the new notation w = —F'. Equation (j9]) for w takes the simple form 
describing an unharmonic oscillator: 



{dl - A)w + U' {w) = 0. 



Potential U{w) is equal to: 
U{w] 



-iT- 2XRo) w + — 
3 ^ 3 



2X 



2nu 



W 



l+2nv 



1 + 2niJ 



:i2) 



TTj) U'{w) = dU/dw. It is useful to remember 
Pm/Pc ^ 1 and hence w -C 1. Thus the first 



where u = l/{2n + 1), g = 2?2A, and in eq. 
that T ^ Rq. Their ratio is about T/Rq ~ 
term in square brackets in eq. (fT2il dominates. Potential U would depend upon time, if 
the mass density of the object under scrutiny changes with time, T = T{t). 
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If only the dominant terms are retained in equations (ITT]) . ( I12p and if the space deriva- 
tives are neglected, equation (ITT]) simplifies to: 

- + r/3-^^ = 0- (13) 
It is convenient to introduce dimensionless quantities: 

t = 7r, w = /3z , (14) 

where P and 7 are so chosen that the equation for z becomes very simple: 

z" - z-" + {1 + kt) = 0. (15) 

Here prime means differentiation with respect to r and the trace of the energy-momentum 
tensor of matter is parametrized as: 

T(t) = ro(l + kt) . (16) 

Constants 7 and /3 are equal to 

^ (-fi„) \ tJ • '■"> 

/ p \ 2n+l 

/3 = 7^To/3 = g (^-^J . (18) 

Thus /9 is a small dimensionless number, and 7 has dimension of time. It is essential that 
7, which determines characteristic time scale, may be much shorter than the universe age, 
t\j, due to the small factor [R^jT^^^^ . Assuming that 3g ~ 1 and i?o ~ lAc/' fi^'i 
n = 2 and pm = 10~^^g/cm^: 7 ^ 400 sec. It would be much smaller for larger n or pm- 
For example if n = 3 and the same pm we find 7 = 0.004 sec. 

In the case of constant T {k = 0) or very slowly varying T [k <^ 1) the solution of eq. 
([T^ is evident. If the initial values z{0) and ^'(0) are sufficiently small, z{t) oscillates near 
the minimum of the potential, which is situated at 

Zmin = (1 + f^r)-'/" . (19) 

If by some reason the magnitude of z{0) takes a sufficiently large value, z > {1 — uY^^ , 
such that potential 

U{z) =z- z^-'l{\ - v) . (20) 

becomes positive, evidently at some stage z{j) would overjump potential U[z^ which is 
equal to zero at 2; = ("the waves are cresting over"). In other words, z(t^ would reach 
zero, which corresponds to infinite i?, and so the singularity can be reached in finite time. 
Analogous situation can be realized if the initial velocity, -2'(0), is sufficiently large. 

The singularity can be also reached in finite time even if z was initially situated at 
the minimum of the potential and the initial velocity was zero. It would take place if 
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Figure 1: Ratio z{t) j ZmmiT) (left) and functions z{t) and Zminij) (right) for 
n = 2, K = 0.01, pjp, = 10^ 



K is positive, i.e. the energy density rises with time. The motion of Zmin to zero and 
simultaneous diminishing of the depth of the potential well make it easier for z{t) to reach 
zero. On the other hand, it is not evident that for decreasing energy density z{t) initially 
resting at the minimum of U{z) could reach zero, most probably it would not. 

We have solved equation ffTSl) numerically and found that indeed the singularity, R — )■ 
oo, is reached in finite time for rising T(r) under quite general conditions. The solution 
for n = 2, K = 0.01, and Pm/ Pc = 10^ is presented in Fig. 1, where the ratio z{t) / Zmin{T) 
(left) and functions z{t) and ^mm(T) separately (right) are depicted. The initial conditions 
are taken as z{0) = 1 and ^'(0) = 0. 

In Figs. 2 and 3 the same quantities are presented for n = 3 and n = 4 respectively. It 
is clearly seen that z{t) reaches zero after a finite number of oscillations around Zmini'^)- 
When Zmin{T) shifts to smaller values, function z{t) initially remains behind but when the 
displacement from the equilibrium point becomes large enough, z{t) started to run after it 
with an increasing speed, then overtakes the position of the minimum, and oscillates back. 
After a few oscillations the retarded position of z{t) happens to be above the point zo{t), 
where the potential is zero. It is essential that the position of this point moves to smaller 
values with rising T(r). Because of that it is easier to overjump the potential at 2; = 0. It 
is intriguing that the magnitude of the ratio z/zmm is approximately equal to 3 in the last 
maximum before the singular point z = is reached. We have not found an explanation 
for that. 

In terms of physical time, t, the evolution of the energy density can be presented as 
T{t) = To(l + t/tch), where tch is the characteristic time of the variation. Coefficient k in 
equation ([TBI) is expressed through t^h as: 

K = ^/t,h . (21) 

Thus the presented in Fig. 1 case of k = 0.01, = 2, and pm/ Pc = 10^ corresponds to 
tc/i = 4 ■ 10^ sec, while n = 3 (Fig. 2) corresponds to tch = 0.4 sec. The characteristic time 
of the density variation can be estimated as tch ~ d/v, where d is the size of the system 
and V is the velocity of the constituent particles in the process of the collapse of the cloud 
or in the collision of the clouds. In the first case the velocity would be quite low and the 
characteristic time is expected to be close to the Newton free-fall time but in the case of 
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colliding clouds the velocities are typically galactic ones, about 300 km/sec. The velocity 
may be even larger at the collision of the supenova ejecta with galactic or intergalactic 
clouds. 

It is more informative to act another way around, namely to estimate k knowing size, 
(i, of the object with changing mass density or sizes of the colliding objects: 



For n > 2 and astronomically large clouds one should expect k <^ 1. For very small k our 
numerical calculations with quickly oscillating functions are not reliable but it seems nat- 
ural to expect that the system would reach singularity according to the analysis presented 
above. 

As it is seen from the numerical calculations, the singularity is reached when t ~ tch- 
This is much shorter than the cosmological time for clouds of denser matter in galaxies or 
a collapsing cloud forming a star or another denser body. 

In the analysis of eq. (ITT]) the spatial derivates have been neglected. At first sight, the 
account of these terms could inhibit formation of singularity, as e.g. happens in the process 
of structure formation due to gravitational (Jeans) instability. However the situation is 
opposite here and the inhomogeneities stimulate singularity formation. Indeed, the effect 
of inhomogeneities can be described by an appearance of the term w/d"^ in eq. ffTTj) with 
positive coefficient. Such a term is equivalent to an addition of an extra attractive force 
pushing w or z to zero, i.e. to — j- oo. 

There are several possible cases when the conditions leading to singularity can be 
realized: collapse of gas cloud leading finally to star formation, collision of two gas clouds 
in a galaxy, stellar ejecta colliding with interstellar or intergalactic matter, and many 
others. From the calculational point of view such processes can be either adiabatic, when 
the mass density changes slowly (this is the case analyzed above) or fast, when the mass 
density changes instantly in an explosive way. Seemingly the latter would result in a faster 
approach to singularity. This case will be analyzed elsewhere. 

If R becomes large, the approximation of fiat space-time would be invalid and the 
derivatives in the equations of motion should be changed into covariant ones. We have not 
analyzed if in this case the approach to singularity is terminated. However, even if it is 
terminated, it takes place at high curvatures when gravity becomes strongly different from 
the Newtonian one. 

Another possible way to avoid singularity is to introduce i?^-terms into the gravitational 
action: 



where m is a constant parameter with dimension of mass. 

With such an extra term in the gravitational action it becomes impossible to express 
analytically R through F'{R). So one needs to work with the equation of motion for R. 
In the homogeneous case and in the limit of large ratio R/Ro equation iQ is modified as 



SF{R) 



R^l^m^ , 



(23) 



1 - 



i? - (2n + 2) 



te R^'^+^R + T) 



0. 



(24) 



6An(2n + l)R\ 




R Q\n{2n + l)Rl 
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To analyze eq. let us introduce, as is done above, dimensionless curvature and 

time: 



R 
J-o 



rp2n+2 -\ 1/2 



Q\n{2n + 



(25) 



Correspondingly eq. ( 12^ is transformed into: 



(1 + gy'-^') y" - 2(n + 1) + y^^^' [y - (1 + k.t,)] = , (26) 
where prime means derivative with respect to ri and 

rp2n+2 

^""6An(2n + l)m2i?2«+i (2'^) 

For very large m, or small g, when the second term in the coefficient of the second deriva- 
tives in eqs. and (12^ can be neglected, the numerical solution demonstrates that R 
would reach infinity in finite time in accordance with the results presented above. Nonzero 
g would terminate the unbounded rise of R. To avoid too large deviation of R from the 
usual gravity coefficient g should be larger than or of the order of unity. Notice that the 
factor (1 + gy'^"'^'^) is always non-zero because (7 > 0. 

Keeping in mind the bound on m > 10~^'^ eV, which follows from the laboratory tests 
of gravity [13], we find n > 6, demanding that the gravity of objects with p ~ 10~^^g/cm^ 
is not noticeably distorted. In ref. [2] a stronger bound is presented, m 3> 10^ GeV. If this 
is the case, then n > 9. A natural value is m ~ mp/ and correspondingly n > 12. For 
smaller values of Tq the bounds on n are noticeably stronger. 

As follows from eq. (12^ . the frequency of small oscillations of y around yo = 1 + kiTi 
in dimensionless time ti is 



i^yl^^i 
l + rf+2 - g 



It means that in physical time the frequency would be 



1 fToV^'^jr^^ 



In particular, for n = 5 and for a galactic gas cloud with Tq/Rq = 10^, the oscillation 
frequency would be 10^^ Hz ^ 10^^ eV. Higher density objects e.g. those with p = I g/cm^ 
would oscillate with much higher frequency, saturating bound ( l29l) . i.e. ~ m. All kind 
of particles with masses smaller than m might be created by such oscillating field. 

On the other hand, as we have seen above, for denser objects the variation of T in 
terms of r or ti is very slow because of very small k. As a result the amplitude of the 
oscillations around the equilibrium point would be also small and possibly such oscillations 
are of no danger from the observational point of view. Still it is possible that there might 
be intermediate cases when the oscillations would lead to observable phenomena. 
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Thus we have shown that the impact of the considered above versions of modified 
gravity on the systems with time dependent mass density in the contemporary universe 
could be catastrophic, leading to the singularity i? — t- oo during finite time in the future. 
This time is typically much shorter than the cosmological one. The problem can be fixed 
by the i?^-term if the power n is sufficiently large, n > 6 (or maybe n > 9). So either 
the versions of the theory with large n or theories with another form of F{R) should 
be considered. A more exciting possibility is that the explosive phenomena predicted by 
modified gravity are observed in the sky. 
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